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We suggest a mechanism whereby the three generations of quarks and leptons correspond to 
surface modes in a five-dimensional theory. These modes arise from a nonlinear fermion dispersion 
relation in the extra dimension, much in the same manner as fermion surface modes in a topological 
insulator or lattice implementation of domain wall fermions. We also show that the topological 
properties can persist in a deconstructed version of the model in four dimensions. 
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I. INTRODUCTION 

It remains a mystery why there are three particle gen- 
erations in the standard model and why they have the ob- 
served pattern of masses and mixing angles, despite many 
attempts at explanation, experimental evidence for flavor 
physics beyond the standard model being limited to neu- 
trino masses. The smallness of neutrino masses and the 
absence of flavor-changing neutral currents and electric 
dipole moments all suggest that the origin of flavor lies 
enigmatically at very short distance. A well explored 
theoretical program is to assume the existence of three 
generations and guess at textures for the mass matrices 
which can be justified by a hierarchy of flavor symme- 
try breaking — an ambiguous exercise given the lack of 
experimental flavor probes in the right-handed fermion 
sector. After pioneering work involving abelian flavor 
symmetries [T], numerous models were also introduced 
with non-abelian flavor symmetries possessing three di- 
mensional irreducible multiplets to justify the existence 
of three generations of quarks and leptons. While this 
general approach can boast of qualitative successes, no 
models have emerged that are particularly compelling. 

Composite and extra dimension models are a natu- 
ral place to look for an explanation for flavor and the 
number of generations: both generically contain towers 
of states, and one can arrange that only three genera- 
tions are light. Composite models must typically rely 
on gauge dynamics to explain the origin of three genera- 
tions, as in Ref. [5], while extra dimension models often 
rely on the Dirac equation having three zero modes in 
certain background flelds of nontrivial topology Ref. [3l- 
[5]. In composite models, the Yukawa matrices of the 
standard model are due to complex interactions between 
constituents, and at best their texture can be predicted; 
in extra dimension models, the Yukawa matrices can be 
computed from wave function overlap integrals in the 
transverse space (see, for example, |51fH]). To the ex- 



tent a gauge/gravity duality pertains, it is possible that 
these two very dissimilar descriptions could be related. 

In this Letter we consider an interesting phenomenon 
observed with lattice domain wall fermions, where the 
number of massless fermions bound to the surface of a 
semi-infinite fifth dimension depended discontinuously on 
the fermion dispersion relation, and hence on coupling 
constants in the action [HI [16]. It was subsequently 
shown that the number of light "families" could be un- 
derstood as a topological property of the five-dimensional 
(5D) fermion dispersion relation in momentum space 
|17j . That is because the number of 4D massless sur- 
face modes is directly related to the quantized coefficient 
of the Chern-Simons operator obtained by integrating 
out the heavy bulk fermions, following the analysis in 
Ref. [IH] ; and this coefficient is obtained from a one- loop 
Feynman diagram which computes a momentum-space 
winding number associated with the fermion propagator. 
This phenomenon was first discussed in the classification 
of fermion modes in liquid helium [19j and is the same 
phenomenon that defines topological insulators [20^221 . 
We consider here that the replication of quark and lepton 
families we observe in the standard model arise in such a 
manner (see [53] for related speculations); an attractive 
feature of the mechanism is that while the number of light 
families is determined topologically, their transverse wave 
functions in the extra dimension are all different and dy- 
namically determined, allowing interesting mass mixing 
without overly restrictive family symmetries. As topol- 
ogy in momentum space depends on the large momentum 
behavior of the fermion dispersion relation, such models 
are forced to confront UV physics, and cannot simply 
rely on an effective field theory description. Therefore 
after describing the general mechanism and providing a 
phenomenological toy example, we look at various UV 
completions that can give rise to a well-defined low en- 
ergy theory. 



II. MULTIPLE ZERO MODES 
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We start by considering fermions in 5D, with an in- 
verse Euclidian propagator which respects 4D Lorentz 
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invariance 

lG-\p^,P5) = iZ.ipYl^ + »Z5(P5)7' - 5](p,P5) , (1) 

corresponding to a plane wave exp(ipaa;°), where 
Pa = {Pi_i,P5\ is the 5-monientum and m is a Dirac spinor. 
We assume Hermitean gamma matrices with and Z5 
being real, odd functions of momentum, and E a real, 
even function, so that corresponds to Hermitean 

derivative interactions in the fermion action. will 
typically have one or more zeros for real p^ and some 
complex value for ps, and a u spinor which is an eigen- 
state of 75. For the generic case ^[ps] ^ 0, this pole in 
G implies a wave function growing exponentially in one 
of the X5 directions; because of the symmetry of Z and 
S under — >■ — P5, if there is chiral solution to = 
with one sign for ^[ps], there is also a solution with the 
opposite chirality and the opposite sign for ^[ps]. Such 
solutions are not normalizable; however, if translation 
invariance in is destroyed, either by a boundary or 
explicit localized 2:5 dependence in the action, it is pos- 
sible that some chiral solutions are normalizable and re- 
tained in the Hilbert space — confined to the boundary 
or defect — while modes of the opposite chirality remain 
non-normalizable and are discarded. The result is an ef- 
fective 4D theory at low energy with chiral fermions; this 
is just a generalization of the domain wall fermion zero 
mode discovered by Jackiw and Rebbi [21] . If there are 
multiple solutions to G^^ = 0, then the 4D theory will 
have multiple families, distinguished by different trans- 
verse wave functions in the extra dimension. While these 
wave functions will depend continuously on parameters 
in the 5D Lagrangian, the number of families can only 
depend discontinuously on those parameters; this makes 
the number of families look like a topological number, 
as indeed it is [TT]. For example, for 5D lattice domain 

wall fermions one finds E = [to -I- rJ2a=ii'^'^^Pa ~ 1)]' 
Za = sinpa, r being a free parameter. It was shown in 
Rcf. 16 that when the extra dimension is made semi- 
infinite, chiral zero modes exist at the boundary, where 
the number of families jumps through the binomial coef- 
ficients {1, 4, 6, 4, 1} as the ratio (m/r) is tuned through 
multiples of two, alternating in chirality with each jump. 
We suggest here that the three families observed in na- 
ture might arise in such a manner. 

A. A toy model 

For a toy model with three generations we start with 
the simple and unjustified assumption that the dispersion 
relation Eq. ([T]) is given by 

iG-^ = p^-f^ + zp5(l + cipl)-f5 - to(1 + C2pI) (2) 

where m, ci and C2 are real and chosen so that for p^ — 0, 
G~^ has three roots, all of a given chirality, occurring at 
P5 — IK with Ki = a, K2 — (b + ic), — [b — ic), where 
a, 6, c are real, positive numbers. These three roots cor- 
respond to transverse wave functions for the zero modes 



of the form (jiiiu) = G^T>{—i^iy)- For now we consider 
the extra dimension to be semi-infinite, and we ignore 
gravity. 

Expanding the 5D fields in this non-orthonormal 
set of transverse zero mode wave functions yields the 
low energy 4D theory with off-diagonal kinetic terms, 
Zijipln^a'^ipj, where ipii^) is the chiral spinor associ- 
ated with the i*^ zero mode, i,j — 1,2,3, and the wave 
function mixing matrix Z is given by the overlap of trans- 
verse wave functions, 

^ I dy 4,* {y) 4>, (y) = . (3) 

Jo '^i + 

4D gauge invariance requires there to be a y-independent 
mode for the W bosons, and so weak currents and Z will 
both be diagonal in the fiavor basis. This is more appar- 
ent if one discretizes the extra dimension and considers 
it as a flavor index, so y-independence of the W is equiv- 
alent to the statement that the W couples to each flavor 
in the UV theory with the same coupling constant. 

The natural starting point for construction of a low 
energy description of the world is to assume one Dirac 
field in the 5D theory for each Weyl field of the standard 
model: Q,U,D,L,E with SU{3) x SU{2) x U{1) quan- 
tum numbers (3, 2)i/6, (3, 1)2/3, (3, l)-i/3, (1, 2)_i/2 and 
(1, l)_i respectively. We assume that Q and L only has 
left-handed chiral zero modes qi, £i, while U, D and E 
have only right-handed zero modes, Ui, di, e^. These 
fields can then couple to the Higgs in the 5D theory as 
UHQ, etc. (Alternatively one could replace U, D, E 
by conjugate fields V^, D'^ and E'^ with Higgs couplings 
WHC^Q, etc, where C5 is the 5D charge conjugation 
matrix, and assume a 5D dispersion relation where all 
zero modes are left-handed.) 

A problem with the above implementation is that, with 
a single Yukawa coupling for each 5D field, it is not pos- 
sible to include weak CP violation. We therefore make 
the Higgs sector the origin of CP violation in the UV 
theory, which requires introducing additional scalars. We 
assume a two-doublet model with a relative CP-violating 
phase in their vacuum expectation values (vevs) arising 
from explicit CP-violation in the Higgs potential. To 
avoid large flavor changing neutral currents, the theory 
is given a discrete symmetry (softly broken) to ensure 
that Hu and Hj^ couple solely to up-type and down-type 
quarks respectively |25j . Furthermore we must assume 
that these are bulk fields and that Higgs vacuum expec- 
tation values are y-dependent; this is the same model and 
mechanism commonly used to introduce CP violating 
bubble walls in theories of electroweak baryogenesis [25] . 
The Yukawa interactions of the quarks in the 5D theory 
are given by {yjjUHuQ + yoDHdQ + h.c). In principle 
there could be higher derivative Yukawa-like interactions, 
but we assume they are zero at the UV scale. 

As with the fermion zero modes, we assume an expo- 
nential form for the 5D profile of the Higgs vevs. By 
means of a hypercharge gauge transformation we can 
make Hd vev to be real and put all the CP violat- 
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ing phase into the {Hd) = ue"^ sin/3, {Hy) = 

yg-(hr+ihi)y where we have chosen the scale of the 

coordinate y so that the exponent in the profile of {Hj) 
equals one, while the real parameters and hi charac- 
terize the profile of Integrating over the coordinate 
y then gives rise to conventional 4D mass matrices, such 
as [Mu]ab oc {yuv cos (i/^){KQ^a + i^u,b + K + ihi)~^. 

In our model we choose Higgs couplings to the leptons 
of the form 
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yELHuE + -{LHaY C^{LHd) + h.c 



(4) 



where H — (J2H*; A has dimensions of mass and controls 
the size of the resulting Majorana neutrino masses. After 
integrating over the coordinate y the Yukawa interactions 
give rise to mass matrices which are simple functions of 
the phenomenological n parameters; from these and the 
wave function matrices Eq. ^ it is straightforward com- 
putation to determine all the masses and mixing angles. 

The point of this toy model is not to present a full 
theory of 5D physics, but only want to see whether 
a model based on this topological insulator mechanism 
could agree with experimental data on flavor physics to 
high accuracy. The model described has 21 real parame- 
ters: a, 6 and c for the Q, U , D, L and E fields, hr and hi 
for the three Yukawa couplings yu, yjj, yE, and the 

scale A characterizing the neutrino coupling to the Higgs. 
These are fit to 18 data: the nine quark and charged 
lepton masses, three CKM angles and one phase, two 
neutrino Am? values, and three neutrino mixing angles. 
Ignoring uncertainties in the data, in general one would 
expect some number of disconnected three-dimensional 
manifolds in parameter space where the model agrees 
with all the data — that number possibly being zero; 
lower dimensional spaces of solutions would generally re- 
quire fine tuning. In searching for solutions, we ignore 
radiative corrections in our model (such as running of the 
Yukawa couplings and masses, or radiative generation of 
higher derivative operators), and fit the parameters to 
data currently available from the Particle Data Group 
[27j . augmented with recent evidence for nonzero ^13 in 
the neutrino sector [28l |29] . 

One might think that this model, with more parame- 
ters than data, could not be predictive. However, when 
we numerically map out the manifold of solutions con- 
sistent with the data (including experimental and theo- 
retical errors quoted in the PDG) we find (i) we always 
get a normal, nondegenerate hierarchy, with 777.3 in the 
narrow range 0.048 eV < 7^3 < 0.051 eV; (ii) solutions 
do not favor maximal mixing for 1^2 ~ '^s; (iii) we find 
in the narrow range -0.023 < Ju < -0.014. What 
has happened is just that the three-dimensional manifold 
of predictions form this constrained model maps onto a 
narrow range of physical properties, so that it is in fact 
somewhat predictive. Nevertheless, one would like to un- 
derstand how to construct a model which is both more 
predictive and is well defined in the UV. 



III. UV COMPLETION: LITTLE FLAVOR 

A general method for constructing a UV completion for 
extra dimension models is to discretize the extra dimen- 
sions while keeping the 4D world continuous (deconstruc- 
tion) [Sni EI] ■ In such an approach, the extra coordinate 
for bulk fermions essentially becomes a flavor index. A 
4D deconstructed version of a free theory with dispersion 
relation similar to Eq. ^ is readily obtained by discretiz- 
ing an infinite extra dimension with defect at site 77 = 0. 
In this case we have an infinite number of flavors of 4D 
fermions with the mass matrix £,„ = [ipi^Mipj^ + h.c), 
where M is an infinite matrix representing the fifth di- 
mension differential operator. For example: 



(Mv) 



icVn-3 + bVn_2 + aVn^l + Vn) 77 < 

[(c7;_3 + bv^2 + + vq) 

+{vo + avi + bv2 + CV3)] 77 = 

^ {Vn + av,i+i + bVn+2 + C^n+s) »^ > 

^1 



(5) 



has the solution Mv — with Vn = x'"', where x is any 
of the three roots of the equation 1 + ax + bx^ + cx^ = 0. 
The vector is normalizable if < 1, and three normaliz- 
able zero modes may be found over a range of parameters 
a, 6, c. The topological nature of the underlying theory is 
manifested by the fact that the the number of normaliz- 
able zero mode solutions does not change with small local 
excursions of M away from the above form. When such 
a system is gauged, however, the gauge fields couple to 
an infinite number of flavors in such a theory, and so the 
theory has a Landau pole and is ill-defined. To cure this, 
one must work with a finite discretized extra dimension. 

With finite continuous extra dimensions, both expo- 
nentially growing and falling zero mode solutions are nor- 
malizable and to obtain a chiral theory one must per- 
form an orbifold projection [321 133j . A well-known ex- 
ample is to compactify the extra dimension as a circle 
parametrized hy 9 € [— tt, tt) with the action respecting a 
Z2 symmetry 'ip{9) — "f^^(—9). This symmetry requires 
mass terms to be odd in 9, and so domain wall defects 
exist at the fixed points of the reflection, 9 — and 9 — tt 
where the zero modes of opposite chirality are located, 
and projecting out modes which are either even or odd 
under this Z2 will result in a chiral theory. A true UV 
completion of our topological insulator model for families 
is possible if a discretized version of this orbifold projec- 
tion exists and is compatible with the survival of multiple 
chiral zero modes. In this case M is a finite dimensional 
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FIG. 1: A discretized Z2 orbifold with three zero modes; TZ 
reflects about the horizontal axis. 
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matrix, and we assume the action respects a Zi symme- 
try ^/j — j^TZip, where TZ is the "reflection" operator in 
flavor space with TZ^ = 1. It is possible to prove an index 
theorem 

(A/r - AA-) - -(A/-+ - A/-+) = Tr7^ , (6) 

where Af^ are the number of left-handed zero modes of M 
with Z2 charge ±1, and Af^ is the same for right-handed 
zero modes. Evidently we would like |Tr7?,| = 3 to obtain 
three chiral families, ruling out discretization of the Z2 
orbifold of a circle, and leading us to consider instead 
the discretization of an extra dimension consisting of two 
circles sharing a point. A simple example is the 7-site 
lattice shown in Fig. [T] where we take Cm to be 

il_(i?2 - R7) + L2_{R3 ~Rl)+ L3_(i?4 - R2) 

+Lz{R'j — Rq) + Li{R^ — R3) + Lq{R3 — R5) 
+L^iRi-R3) + h.c. (7) 

where Li and Ri are tpL and 'ipR respectively at site i. 
This has the simple interpretation of hopping terms for 
fermions around the circles, with fermions at the shared 
point (site 3) being able to hop onto either circle. In 
this example, TZ reflects sites about the horizontal axis 
in Fig. [1] and the action is invariant under "0 75 ■ 



The theory has three massless Dirac fermions, and if we 
project out all states with negative Z2 charge, then we 
are left with three left-handed zero modes. The impor- 
tant point of this model is not that M has three zero 
mode solutions, but that the zero modes persist even if 
M is perturbed in any random way which respects the Z2 
symmetry — although the eigenvectors ( "transverse wave 
functions" ) will be altered under this perturbation — ex- 
actly as one would expect from the topological origin of 
this model. 

This mechanism in general — and the deconstructed 
model in particular — naturally suggests a number of in- 
teresting directions to pursue, most interestingly whether 
it imposes inescapable constraints on Higgs and CP 
physics that might be probed by the LHC, as well as 
observable flavor changing neutral currents or lepton fla- 
vor violation, which would be characterized by the scale 
of the extra dimension [M] . 
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